
Collaborative Conference on Crystal Growth (3CG) 

November 2014 


Crystal Growth of Ternary Compound 
Semiconductors in Low Gravity Environment 


Ching-Hua Su 

NASA/Marshall Space Flight Center 


1 


Crystal Growth of Ternary Compound Semiconductors 


A low gravity material science experiment has been prepared 
to be performed in the Low Gradient Furnace (LGF) in the 
Material Science Research Rack (MSRR) on International 
Space Station (ISS). 


There are two sections of the flight experiment: 

(I) Investigation toward crystal growth by physical vapor 
transport (PVT): the growth of ZnSe and related ternary 
compounds, such as ZnSeS, ZnCdSe, and ZnSeTe, 

(II) Investigation on the melt growth of CdTe and CdZnTe by 
directional solidification. 



Crystal Growth Activities at NASA/MSFC 


Technological significance: 

Growth (melt and vapor) and characterization of II-VI 
and IV- VI compounds semiconducting materials, such as 
HgCdTe, HgZnTe (for IR detectors), CdS and ZnO (for 
UV detector), ZnSe, ZnSeTe (for green /blue laser), CdTe 
and CdZnTe (for x-ray, gamma ray detectors), PbTe, 
PbTeSe, and PbSnTe (thermoelectrics). 
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ray telescope 
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II- VI semiconducting compounds grown at MSFC 


Compounds 

HgTe HgCdTe CdTe CdZnTe ZnTe CdS ZnSe ZnS 

Melting 
points (°C) 

670 700 1092 1130 1292 1397 1526 1718 

Melt growth 



PVT growth 
temperature (°C) 


850 


1000 985 1120 1150 


Scientific Objectives 


1 . To establish the relative contributions of gravity-driven fluid 
flows, both in liquid and vapor, to (1) the non-uniform 
incorporation of impurities and defects and (2) the deviation 
from stoichiometry and (3) the compositional variation 
observed in the grown crystals. 

2. To assess the self-induced strain developed during processing 
at elevated temperatures and retained on cooling caused by the 
weight of the crystals. 

3 . The relation between fluid phase processes and the generation 
of defects in a grown crystal is an outstanding problem in 
materials growth. Studies in microgravity will be compared 
with modeling and will lead to a greater understanding of the 
processes involved. 
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(3») X 


Phase Diagram of Zn-Se 



x 


Se 


Partial pressures along the three-phase curve for ZnSe(s) 


At 1 156 °C, P Zn varies from 8.4 to 1 .3x1 (H atm, P St , , varies from 20 
to 5x1 O' 9 atm and a = P Zn / P Se2 varies from 1.7 xIO 9 to 6.5x1 O' 6 




1000/T 


Mass Flux (mole /cm 2 - sec) 


One-Dimensional Diffusion model of PVT 





Calculated mass flux of 
ZnSe as a function of 
AT and different values 
of a(L). The source 
temperature was 
1080°C. Solid lines are 
for a(L) > 2 and dashed 
lines are for a(L) < 2. 


Calculated mass flux of 
ZnSe as a function of 
AT under the same 
conditions except a 
residual gas pressure 
of 0.008 atm is present 
in the system. 


Calculated mass flux of 
ZnSe as a function of 
residual gas pressure 
for source temperature 
at 1080°C and different 
values of a(L). 


Su, et al. (1998); Sha et al. (1995) 



Simultaneous measurements of mass flux and partial pressure 




Simultaneous measurements of mass flux and partial pressure 



Simultaneous measurements of mass flux and partial pressure 


• P Zn measured over the deposit; a (deposit) = 17. 

• The measured P z value was used in the calculation. 



Pressure (atm) 


Calculated partial pressure profiles for ZSTO-3 run 2 



a = 6.4 


Distance (cm) 


Summary of one dimensional diffusion analysis 


• Four experimentally adjustable parameters, the source 
temperature, the deposition temperature, the partial pressure 
ratio over source and the residual gas pressure, determine the 
diffusive mass flux of a PVT system. 

• However, two of these four parameters, the partial pressure 
ratio over source and the residual gas pressure, are more 
critical than the others. As will be shown, these two 
parameters are critically dependent on the proper heat 
treatments of the starting materials for optimum mass flux. 
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Optimum Heat Treatments of Starting Materials 


P Se in equilibrium with ZnSe for various optical cells 


H 2 reduction and vacuum treatments for ZnSe(0)-3 established to be the optimum procedure 

1145 °C 1035 °C 



1000/T (K) 


Su et al. (1998) 
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Crystal Growth by Physical Vapor Transport: 
temperature profile and initial ampoule positions 


• The growth ampoules can be equipped with optical windows to confirm the 
stoichiometry of the starting material before growth. 

• The thermal profile, with a maximum in the middle, was provided by a three-zone 
furnace with an adiabatic zone between central and cold zones. 

• Growth was initiated by translating furnace to the right. Gravity 



Distance 


Su (1995) 
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PVT Growth of ZnSe and Related Ternary Compounds 


1 . Self-seeded growth of ZnSe in vertical (stabilized and 
destabilized) and horizontal configurations 

2. Seeded growth of ZnSe in vertical and horizontal 
configurations 

3. Self-seeded growth of Cr-doped ZnSe in vertical and 
horizontal configurations 

4. Self-seeded growth of ZnSeTe in vertical and horizontal 
configurations 

5. In-situ and real-time optical monitoring of seeded growth 
in a horizontal configuration 
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Flow chart of sample characterization plan 
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Gravity Effects on the Grown Crystals 


Effects were studied by comparing the following 
characteristics of horizontally and vertically grown ZnSe 
crystals in : 

• Grown crystal morphology : contactless growth for the 
horizontal configuration. 

• Surface morphology of the grown crystals was examined by 
SEM and AFM. (growth was terminated by stopping furnace 
translation, lowering the source temperature by 1 0 °C and then 
cooling the thermal profile at the same rate) 

• Segregation and distribution of defects and impurities in the 
grown crystals was determined by photoluminescence, SIMS 
and precision density measurements. 
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Gravity Effects on the Grown Crystals 


Morphology of the as-grown crystals: 

I. Self-seeded ZnSe: Crystals grown in the horizontal configuration grew away 
from the ampoule wall with large (110) facets tend to align parallel to the 
gravitational direction. Crystals grown in the vertical configuration grew in 
contact with the wall to the full diameter. 

II. Seeded ZnSe: the as-grown seeded crystals for the horizontal and vertical 
configurations showed similar characteristics in the morphology as described 
above for the self-seeded growth. 

As-grown surface morphology: 

I. As-grown surface of horizontally grown crystals was dominated by (110) 
terraces and steps (identified to be (221) in one case). 

II. As-grown surface of the vertically grown 

(a) Crystals showed granular structure with nanotubes (200nm OD, 75nm ID, 
25nm in height for one case on ZnSe) on the top. 

(b) Some crystals showed a network of high plateau with each island 30 - 
70mm in diameter and 3.5mm in height. Numerous nuclei were observed with 
diameter 20 - 50nm and height of 1 - 7nm on top of these islands. 
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Morphologies of Self-seeded Vertically 
Grown ZnSe Crystals 


ZnSe-25: vertically stabilized ZnSe-3 1 : vertically destabilized 

ampoule ID: 15mm _ 



t 
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Morphologies of Self-seeded Horizontally Grown ZnSe Crystals 






Gravity Effects on As-grown crystal morphology 


• For furnace translation rates higher than the mass flux: 

- In the horizontal configuration the crystal maintained the growth surface 
morphological stability by (1) self-adjusting the degree of supersaturation to 
increase the mass flux or/and (2) by reducing the cross section area of the 
grown crystal. 

- In the vertical configuration the crystal growth surface became 
morphologically unstable with voids and pipes embedded. 


ZnSe-47H 7.5 mm/day ZnSe-35V 11.4 mm/day 



Morphologies of Seeded Vertically Grown ZnSe Crystals 



ZnSe(S)-9V ZnSe(S)-12V 



23 




Morphologies of the Seeded Horizontally Grown ZnSe crystals 


x 
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ZnSe(S)-13H 


ZnSe(S)-8H 


Gravity Effects on Impurities and Defects Distribution 


I. Results from SIMS mappings : 

For the horizontally grown self-seeded ZnSe crystal [Si] and [Fe] showed clear 
segregation toward the bottom on the wafer cut axially along the growth 
axis. 

For the vertically grown seeded ZnSe crystal [Si] and [Cu] showed segregation 
toward the edge of the wafer cut perpendicular to the growth axis . 

II. Mappings of PL near band edge intensity ratios indicated: 

(1) All the horizontally-grown crystals showed the following trends in the 
radial and axial segregation of [Al] and [ V Zn ] : 

[Al] segregates radially toward the top and axially toward the first grown 
region and [V Zn ] segregates radially toward the bottom and axially toward the 
last grown region. 

(2) The as-grown surface of the seeded vertically stabilized grown crystal 
showed [Al], [Li or/and Na] and [V Zn ] segregate radially toward the center. 

(3) The as-grown surface of the self-seeded vertically destabilized grown 
crystal showed [Al] and [V Zn ] segregate radially without an apparent pattern. 


Su et al (1999) 
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x (mm) 


Impurities Distribution by SIMS (horizontally grown) 



13 - 



11 - 



cone. Si (ppb) 
2.04312 
1.86626 
1.00937 
1.0125 
1.73562 
1.65075 
1.50107 
1.505 
1.42012 
1.35125 
1.27437 
1.1975 
1 . 1 2062 
1.04375 
0.966075 


y (mm) 
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13 


12 
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Concentration of Fe In ZnSe-40 


all paints 



cone. Fe (ppb) 

2556.75 

2390.5 


2050 


1725.5 

1559.25 


1226.75 


I 


720 


395.5 


J i i i i I i i i i I i i i i I i i i i I i i i i I l_ 

5 6 7 8 9 10 

y (mm) 


Concentration of Cu In ZnSe-40 



y (mm) 
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Impurities Distribution by SIMS (vertically grown) 


141 — 


ncentratlon of SI In ZnSe- 16 v 



cone. Si (ppb) 

5.56937 

6.23675 

5.60612 

5.3775 

4.94567 

4.51525 

4.06562 

3.555 

3.22437 

2.79375 

2.36312 

1.9325 

1.50167 

1.07125 

0.540625 


gravity 


all paints 


I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

6 7 8 9 10 11 12 13 14 

x (mm) 
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gravity and 
growth direction 


141 - 


n ZnSe-IGv 


cone, of Fe (ppb) 

H ■ 



12526 

11665.5 
1 1203 

10540.5 
9676 

9215.5 
6553 

7690.5 
7226 

5565.5 
5903 

5240.5 
4576 

3915.5 


10 11 

x (mm) 


Hr 
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12 


1 1 

'e 

E 10 


9 - 


8 - 


7 - 


A % 

Concentration of Cu In ZnSe- 16 v 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


cone, of Cu (ppb) 
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731.166 
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517.436 
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Summary on the Photoluminescence I 2 and Ij dee P emissions 


l 2 emission: 

• l 2 , the exciton bound to substitutional donor, emission in our ZnSe samples can be 
attributed mainly to Al impurity, with A(I 2 )/A(F X ) = 4.88 corresponding to 1700 ppb, 
atomic, or 7.46x1 0 16 cm -3 . 

• Isshiki et al. (1991) gave the expression between intensity ratio (l 2 )/(l Fx ) an( J N D : 

log 10 (l 2 /l Fx ) = -22.0775 + 1.46268 log 10 N D (cirr 3 ) 

• Therefore, (l 2 /l Fx ) = 82 [A(I 2 )/A(F X )] 


\^deep emission: 

• I/^p is related to exciton bound to V Zn deep acceptor and [V Zn ] is proportional to 
A(l 1 dee P)/A(F x ). 

• The reaction during Zn vapor annealing: 




Zn zn V Zn + Zn(g) 




K = [V Zn ] x P Zn = K 1 A(l 1 dee P)/A(F x ) x P Zn 
The ZnSe samples were annealed at 1104 °C: 

1. A(l 1 dee P)/A(F x ) = 7.52 when sample is in equilibrium with P Zn 
= 6. 1x1 0 -3 atm (a = 6.05) A(l 1 dee P)/A(F x ) x P Zn = 0.0459 atm. 

2. A(l 1 dee P)/A(F x ) = 5.18 when sample is in equilibrium with P Zn 
= 9.0x1 O' 3 atm (a = 19.43) Afl^pyAfFJ x P Zn = 0.0466 atm 



2.77 2 775 2.78 2.785 2.79 2.795 2.8 2.805 

Photon Energy (eV) 
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Distribution of [V Zn ] and [Al] in ZnSe (horizontally grown) 





1 1 I I I II I 1 I I I I I I I I I I I I I I 1 I I I I II I I I I I 11 I I 1 1 1 1 11 1 l,.l I I I I I I 11 
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x distance (mm) 



A(I?)/A(F x ) 
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26.2654 

24.9909 
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growth 
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Distribution of [V Zn ] and [Al] in ZnSe (horizontally grown) 



t 1 1 I L 

0 2 4 6 8 

distance, x 


A(I 2 )/A(F x ) 



] gravity 


11.679 

11.279 

10.880 

10.481 

10.082 

9.6825 

9.2830 

8.8838 

8.4845 

8.0853 



I , , L 

10 12 


I l l l 

0 2 4 6 8 

distance, x 


gravity 


X 

growth 

direction 


A(V)/A(F X ) 

39.678 
37.448 
35.219 
32.990 
30.761 
28.531 
26.302 
24.073 
21.844 
19.614 



| gravity 


I L 

10 12 
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Distribution of [V Zn ] and [Al] in ZnSe (vertically grown) 
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Vertical stabilized 
configuration 


A(I, d )/A(F x ) 
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a ( i ;)/ a ( f x ; 

1.003 
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0.859 

0.811 
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• growth direction 
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Distribution of [V Zn ] and [Al] in ZnSe (vertically grown) 


Vertical destabilized configuration 
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Synchrotron X-ray Transmission Topography of ZnSe wafer 


Thin twin lamellae (TL) and three sets of (111) slip bands (S) originating 
from the lateral twin boundaries. 




(diffraction 

vector) 


x gravity 


ampoule 

wall 


g = (331) 


2mm 


g = (133) 


g = (313) 



Reflection Topograph of ZnSe-43H as-grown Facet 


Reflection SWBXT image on the as-grown facet of ZnSe-43H showed 
areas of twins. The crystalline quality is generally good except that the 
upper region, where the crystal started to grow into full diameter and 
away from the ampoule wall, showed lattice strain. 



1mm 




ZnSe 43H 


In-situ optical monitoring during crystal growth 


• Optical absorption for partial pressure measurements along 
the length of the growth ampoule to measure vapor phase 
transport species distribution 

• Optical interferometry to measure seed thermal etching, 
instantaneous growth rate and the evolution of surface 
topography 


Visual observation of the growth evolution 


Ampoule Design for in-situ Optical Monitoring 
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Measured Partial Pressure as a Function of Time 
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2.8 (±.5) mm 


Interferometry Results from Seed Crystal Surface 

(Room Temperature) 



RMS : 0.171* PHASE MAP 14:53:39 

P-U: 1.262* 01-14-98 


Scale : 

1 2.23 
2.14 
2.04 
1.94 
1.84 
1.75 
1.65 
1.55 
1.46 
1.36 
1.26 
1.16 
1.07 
0.97 


Tilt 
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mother plot Cy2N)7 


H 3.7 (±.5) mm ► 


Fringes from Fabry-Perot 
Interferometer 


Fast/Vai Phase Map of Crystal 
Surface, (from Fringes) 





Interferometry Results from Seed Crystal Surface 

(at 1120°C) 



Fringes from Fabry-Perot 
Interferometer 


Fast/Vai Phase Map of Crystal 
Surface, (from Fringes) 



Synchronized plots of sample temperature and the 

measured thermal expansion 



41 


Expansion (mm) 


°in v 


Thermal Expansion Coefficient of ZnSe from Interferometry 


Expansion AL — n(A,/2) X — 632. 8nm L 0 =4. 5 in in 
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Numerical Modeling of Physical Vapor Transport 


• Two dimensional and three dimensional calculations 

• Finite element technique - Fidap code 

• Thermal and Species induced buoyancy forces 

• Compressible or Boussinesq model 

• Benchmark -2D (H 2 -I 2 system - PVT growth) 

• ZnSe calculations with residual gas 

• Benchmark -3D (Natural convection in a cylinder) 

• 3-D ZnSe calculations with residual gas 

• Ongoing and future work 


Ramachandran et al. (2000) 
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The Physical Model 
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Dimensionless wall temperature distribution - PVT 
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Benchmark 


the H 2 -I 2 System 


• benchmark calculations by Rosenberger et al. (J. Crystal 
Growth 51426 1981; 67 241 1984; H8 49 1992) 

• source temperature Ts=370.5 K 

• crystal temperature Tc=358. 1 K 

• ampoule pressure : 1 00 torr 

• I 2 (M=254) is the deposited species and H 2 (M=2.016) is inert 

• 2-d Cartesian system 

• linear wall temperature 

• quasi-compressible and Boussinesq calculations 

• Peclet number analysis, Pe ~ 1 for diffusive flow 

• growth rate results 
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3-D Computational Grid and Code Validation Results 





* Schirolcy and Rosenberger (1984) 
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Parameters of ZnSe System with Residual Gas 


a = 2.9 

• density, 1 . 2x 1 O' 5 g/ cm 3 

• dynamic viscosity, 4.3xl0' 5 Pa-s 

• kinematic viscosity, 36 cm 2 /s 

• diffusivity, Zn in N 2 , 64.59 cm 2 /s 

• diffusivity, Se 2 in N 2 , 71 .46 cm 2 /s 

• thermal expansion coefficient, 7.1xl0' 4 K 1 

• Prandtl number, 0.439 

• Schmidt number, Sc zn = 0.557 

• Schmidt number, Sc se = 0.503 
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) Velocity Difference Plots - procedure 

• calculations for pure diffusion limited conditions, Og 

• calculations for other g conditions, a conditions, etc. 

• calculate differences in axial (u) and transverse (v) 
velocities at all identical grid locations between previous 


• contour the velocity differences and plot 

• determine g-sensitivity 

• Note: all calculations were using the Boussinesq approx. 


steps; e.g. U(0g) - U(lg) and V(0g)-V(lg) 


conventions: 



y,v t 


source 



crystal 

To ; [T=T(y)] 


0=0; horizontal 

0=90; vertical stabilized 

0=270; vertical destabilized 
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Gravity Effect on Velocity Difference (horizontal case) 


- g-effects with constant crystal temperature , Tc. 

AU~ 0.3 mm/s; AV~ 50 pm/s 

- g-effects with crystal temperature variation, Tc = Tc(y). 
AU~ 0.3 mm/s; AV~ 43.75 pm/s 

, Av velocity (m/s); (1 g Q - 0 g Q ) 

2.5 - 

2 - 


Av (m/s) 

4.376E-05 
9.045E-06 
] -8.312E-06 
1 -4.303E-05 
-7.774E-05 
-1.125E-04 
-1.472E-04 
I -1.819E-04 
I -2.166E-04 



- g-sensitivity (horizontal case) based on max. buoyancy 
driven flow normal to growth direction is 1 0% of crystal 
growth rate (3mm/day or 0.035 pm/s) 

transverse acceleration requirement: ~ lxl O' 4 g 0 
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Gravity Effect on Velocity Difference (vertical case) 


- g-effects with crystal temperature variation, Tc = Tc(y). 
Stabilized orientation: AU~ 23.1 jum/s; AV~ 9.4 pm/ s 

- g-effects with crystal temperature variation, Tc = Tc(y). 
Destabilized orientation: AU~ 18.1 pm/s ; AV~ 9.4 pm/s 

As far as transverse velocity difference is concerned both 
vertically stabilized and destabilized orientations have 
similar effects 

- g-sensitivity (vertical case) based on max. buoyancy 
driven flow normal to growth direction is 1 0% of crystal 
growth rate (3mm/day or 0.035 pm/s) 

longitudinal acceleration requirement: ~ 2.7x1 O' 3 g 0 
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ZnSe with Residual Gas - 3D Calculation Results 


• significant flow observed along the ampoule axis (z-direction) 
indicative of more deposition in the central area than near the 
walls 

• velocity contours in the cross planes(x-y) show appreciable 
variation only near the end walls (source and crystal) 

• species (Zn and Se 2 ) show fairly uniform distributions in the 
cross planes 

• predicted crystal growth rate from 2-D and 3-D calculations 
are in fair agreement 
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Summary of Theoretical Calculation 


• 2D and 3D calculations performed for ZnSe system 

• Residual gas effects considered 

• Calculations show that shear flow velocities of 10 to 50 
microns/s are induced by buoyancy effects (290 to 1400 
times growth rate) 

• g-level requirements established based on time scale analysis 
required transverse g level: < 1.2xl0' 4 g o 

required longitudinal g level < 8.5x1 0‘ 3 g 0 

• It is noted that the Boussinesq model used in the calculations 
tend to underpredict velocities 
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Flight Experiments on International Space Station 


• The flight experiments will be 
conducted in the Low Gradient 
Furnace (LGF) in the Microgravity 
Science Research Rack (MSRR) on 
International Space Station (ISS) 

• Nine different growth runs will be 
performed for ZnSe, Cr-doped 
ZnSe, ZnSeTe, ZnSeS and ZnCdSe 
materials with different growth 
parameters. 

• The flight experiments are 
scheduled to commence in late 
2015 . 
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